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1.0  INTRODUCTION  AND  SUMMARY 


Theoretical  work  performed  under  Air  Force  contract  [1]  has  indicated 

that  some  Ga^  alloy  compositions  that  are  lattice  matched  to 

InP  have  electronic  properties  which  give  them  the  potential  for  use  in 

high  frequency  MESFET  devices  which  will  perform  significantly  better  than 

those  made  with  GaAs.  As  an  example,  the  projected  gain-bandwidth  product 

for  a  MESFET  made  with  GaA  T,In_  ,..PA  ,Asa  ,  is  about  2.3  times  that  for 

0.27  0.73  0.4  0.6 

one  made  with  GaAs  using  the  same  device  dimensions  and  doping  level. 

Although  the  growth  by  liquid  phase  epitaxy  (LPE)  of  Ga^  As^ 

lattice  matched  to  InP  [2-4]  and  to  GaP,  As  [5]  substrates  had  been 

i-y  y 

reported  prior  to  mid  1976,  the  compositions  and  doping  levels  used  were 
appropriate  to  optoelectronic  devices  rather  than  to  MESFET  fabrication. 

The  latter  operation  requires  very  thin  (on  the  order  of  0.2  to  0.4  pm) 
epitaxial  layers  of  high  crystalline  perfection.  T^is  implies  an  excellent 
lattice  match  to  the  substrate,  which  must  be  semi- insulating  in  order  to 
avoid  spurious  conduction  and,  furthermore,  must  itself  be  of  high  crystal¬ 
line  perfection  in  order  to  minimize  the  formation  of  electrically  active 
defects  at  the  substrate-epilayer  interface. 

In  August  1976  the  Research  Triangle  Institute  (RTI)  proposed  to  the 
Air  Force  Office  of  Scientific  Research  a  research  effort  to  synthesize  a 
quaternary  alloy  material  and  use  this  material  in  a  prototype  MESFET  to 
enable  comparison  of  performance  with  similar  devices  made  from  GaAs. 

This  research  program  was  envisioned  as  a  three-year  effort  with  the 


following  major  objectives: 


1.  To  demonstrate  LPE  growth  of  the  desired  Ga^_xInxP^_^ASy  alloy 
composition. 

2.  To  use  compositional  and  structural  analysis  to  guide  the 
growth  process  in  order  to  obtain  epitaxial  material  of  the 
highest  possible  perfection. 

3.  To  perform  electronic  properties  measurements  to  verify  theo¬ 
retical  predictions  of  enhanced  properties  relative  to  GaAs. 

4.  To  identify  and  demonstrate  suitable  ohmic  and  Schottky  barrier 
metallization  materials  necessary  for  MESFET  fabrication. 

5.  To  fabricate  and  demonstrate  prototype  MESFET  design  using  the 
the  best  material  produced  and  to  relate  the  performance  of 
this  device  to  material  properties. 

Prior  theoretical  calculations  [6]  had  indicated  that  the  presence  of 

compositional  disorder  in  the  lattice  arrangement  of  the  Ga.,  In  P,  As 

J-“X  x  i  —y  y 

alloys  could  introduce  a  scattering  mechanism  (so-called  "alloy  scattering") 
which  could  reduce  electron  mobility  and,  therefore,  performance  of  devices 
made  from  disordered  alloys.  One  of  the  correlations  sought  in  the  pro¬ 
posed  research  was  that  between  compositional  disorder  and  electronic 
properties. 

This  proposal  was  accepted  by  AFOSR,  and  the  first  year  funding  started 
March  1,  1977.  The  first  year  efforts  have  been  previously  reported  [7]. 
Growth  efforts  during  the  first  year  were  directed  first  to  compositional 
control  to  obtain  alloys  lattice  matched  to  InP  substrates,  particularly 
those  compositions  near  Ga^  73^0  4As0  6'  Theoretical  calculations  [6] 

of  velocity-field  characteristics  had  indicated  that  a  range  of  alloy  compo¬ 
sitions  lattice  matched  to  InP  provided  transport  properties  that  made  them 
potentially  superior  to  GaAs  in  MESFET  fabrication.  The  capability  for 
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growth  of  good  quality  epitaxial  layers  on  (lll)B  oriented  InP  substrates 
was  demonstrated. 


Although  the  quaternary  layers  were  of  good  crystallographic  quality, 

the  electrical  properties  were  not  as  desired.  Even  in  non-intent ionally 

doped  material,  electron  carrier  concentrations  were  large,  initially  on 
18  —3 

the  order  of  10  cm  .  Improved  melt  bake-out  procedures  and  the  use  of 
higher  purity  source  materials  enabled  the  reduction  of  carrier  concentra¬ 
tions  to  the  range  4  *  10^  to  2  *  lO1^  cm  The  highest  electron 

2 

mobility  value  obtained  was  about  4,000  cm  /V-sec  at  a  carrier  concentra- 
16  -3 

tion  of  3.8  x  10  cm  ,  with  both  measurements  made  at  room  temperature. 
This  value  for  mobility  was  about  half  the  predicted  (theoretical)  value. 
The  use  of  InP  buffer  layers  between  the  InP  substrate  and  the  quaternary 
layer  appeared  to  have  no  effect  on  layer  properties. 

Photoluminescence  measurements  for  the  quaternary  alloy  layers  made 
early  in  the  program  showed  the  presence  of  deep  levels,  but  well  before 
the  end  of  the  first  year's  efforts,  these  deep  levels  had  been  eliminated; 
and  the  emissions  from  the  quaternary  alloys  were  intense  band-to-band 
transition  lines  of  narrow  width,  indicating,  along  with  the  improvement 
in  electrical  properties,  that  the  growth  process  was  yielding  much 
better  quality  layers. 

As  a  means  of  measuring  structural  order  (for  correlation  to  elec¬ 
tronic  properties).  X-ray  diffraction  measurements  were  attempted.  This 
required  removal  of  the  InP  substrate  to  allow  direct  exposure  of  the  epi¬ 
taxial  quaternary  layer  to  the  X-ray  beam.  The  fragility  of  the  thin 
(5  to  10  pm)  exposed  epilayers  created  many  problems  in  sample  preparation. 
This  necessitated  the  development  of  a  selective  etch  that  dissolved  the 
InP  substrate  material  at  a  much  higher  rate  than  the  quaternary  alloy 
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material.  With  this  etch  suitable  samples  for  X-ray  analysis  were  being 
prepared  near  the  end  of  the  first  year's  effort. 

Attempts  to  make  Schottky  barriers  on  the  quaternary  material  were 
not  successful  during  the  first  year.  Although  I-V  characteristics  for 
the  structures  formed  by  deposition  of  gold  on  the  quaternary  were  nonlinear, 
there  did  not  appear  to  be  any  well-defined  barrier. 

During  the  second  year  quaternary  growth  was  carried  out  on  (111) A  and 
(100)  oriented  InP  substrates  as  well  as  (lll)B  for  a  total  of  140  growth 
runs.  Many  of  these  rims  were  required  to  work  out  the  phase  diagram,  par¬ 
ticularly  for  the  (111) A  face,  to  enable  compositional  control  over  the 
quaternary  epilayer.  Lattice  mismatch  between  the  epilayer  and  InP  sub¬ 
strate  was,  for  many  layers,  reduced  to  below  the  limit  of  detection 
(±0.03  percent  in  our  case).  Thickness  control  of  about  ±20  percent  was 
achieved  for  epilayers  of  1  to  5  pm  thickness.  For  submicron  layers,  how¬ 
ever,  growth  rate  variations  from  run  to  run  were  much  larger.  Further 
work  is  needed  to  achieve  precise  control.  Distribution  coefficients  for 
the  n-type  dopant  species — tin,  tellurium,  and  selenium — were  determined. 

Photoluminescence  and  X-ray  diffraction  measurements  indicated  high 
crystalline  quality  and  homogeneity  for  those  layers  grown  after  phase 
diagram  data  enabled  precise  composition  control  to  achieve  lattice  match 
between  the  alloy  epilayers  and  InP  substrates.  The  use  of  (111) A  and 
(100)  oriented  epilayer  layers  enabled  relatively  straightforward  prepara¬ 
tion  of  samples  for  X-ray  analysis  to  ascertain  order.  Low  intensity 
diffraction  peaks,  which  could  possibly  be  interpreted  as  due  to  the  pres¬ 
ence  of  a  tendency  to  long-range  order,  were  detected;  but  the  degree  of 
ordering  was  not  quantified.  In  general,  experimental  difficulties  pre¬ 
vented  the  collection  of  enough  X-ray  diffraction  data  to  enable 
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correlation  of  the  degree  of  lattice  ordering  with  growth  conditions  and 
electronic  properties. 

Conductivity  and  Hall  effect  measurements  indicated  that  even  though 

16  —3 

the  non-intentional  carrier  concentrations  were  reduced  to  about  10  cm  , 

there  was  no  improvement  in  carrier  mobility  over  t  *  it  obtained  during  the 

first  year.  Room  temperature  values  of  mobility  in  doped  and  undoped  . 

2 

layers  were  typically  2,000  to  3,000  cm  /V-sec.  This  was  far  below  theo¬ 
retically  predicted  values.  An  analysis  of  mobility  versus  temperature 
data  indicates  that  the  epilayers  probably  contain  a  large  number  of 
inhomogeneities  that  give  large  components  of  space  charge  scattering  and, 
possibly,  dislocation  scattering. 

No  problems  were  experienced  in  making  ohmic  contacts  to  the  quater¬ 
nary  alloys,  but  efforts  to  obtain  stable  and  reproducible  Schottky 
barriers  were  for  the  most  part  not  successful.  Five  different  metals  were 
investigated:  Au,  Al,  Cr,  Ag,  and  Mg.  Of  these,  only  Au  provided  appre¬ 
ciably  nonlinear  I-V  characteristics.  The  largest  barrier  height  obtained 
was  0.51  volt  and  a  diode  factor  n  =  2.6. 

A  preliminary  design  for  X  and  Ka  band  MESFETs  was  undertaken,  but 
because  of  the  low  mobility  values  and  difficulty  in  obtaining  good 
Schottky  barriers,  actual  fabrication  of  devices  was  not  carried  out. 

There  still  appear  to  be  many  problems  to  be  worked  out  before  micro- 

wave  MESFET  device  quality  Ga,  In  P1  As  material  can  be  routinely 

x^x  x  x-y  y 


achieved. 


2.0  MATERIAL  GROWTH 


2.1  Liquid  Phase  Epitaxy 

For  LPE  growth,  standard  horizontal  slider  boats  made  of  high  purity 
graphite  were  used.  Two  different  LPE  reactors  were  committed  to  this 
program  this  year.  One  was  the  Lindberg  "Heavy-Duty"  furnace  described  in 
the  previous  annual  report  [7],  The  second  furnace  was  constructed  at  RTI. 
It  consists  of  "clam-shell"  type  resistance  heaters  containing  a  60-cm 
long  liquid  sodium-filled  heat  pipe  mounted  on  wheels  which  run  on  a  track 
arrangement  that  allows  the  entire  furnace  to  be  pulled  over  the  quartz 
reactor  tube.  Fast  heat  up  and  cool  down  can  be  realized  this  way.  In 
addition,  the  heat  pipe  allows  a  flat  temperature  zone  approximately  25  cm 
long. 

The  LPE  growth  procedure  used  in  most  experiments  was  the  same  as 
that  described  in  the  previous  annual  report  [7] — i.e.,  the  "transient 
method"  with  ramp  cooling  was  utilized  for  most  of  the  growth  runs.  How¬ 
ever,  for  certain  specific  experiments,  to  be  described  later,  the  "step 
cooling"  approach  was  used.  In  these  cases  the  melt  was  saturated  at  a 
given  temperature,  Tg,  (in  the  range  650-640°C)  and  was  then  subjected  to 
a  rapid  cooling  step  of  ATg  varying  between  2  and  10°C  to  cause  super¬ 
saturation  in  the  melt.  In  our  furnaces  such  a  drop  in  temperature  can  be 
achieved  in  two  to  six  minutes.  Once  the  temperature  had  stabilized,  sub¬ 
strate  wafers  were  inserted  in  the  melt  for  the  desired  length  of  time  for 
LPE  growth  to  occur,  typically  ten  minutes. 

Finally,  for  a  few  growth  runs  with  (100)  oriented  substrates,  the 
"two-phase"  method  [2-4]  was  utilized  in  conjunction  with  step  cooling. 

In  this  case  the  melt  was  an  undersaturated  solution  of  Ga,  In,  and  As 


into  which  an  excess  of  solid  InP  was  introduced  in  order  to  supply  the 
required  amount  of  phosphorous  to  the  melt.  With  this  method,  growth 
occurs  on  the  InP  substrate  as  well  as  on  the  solid  phase  InP  present  in 
the  melt. 

In  order  to  facilitate  analyses  by  X-ray  techniques  to  determine  the 
state  of  ordering  in  these  quaternary  compounds,  it  was  found  desirable 
that  epilayers  should  be  grown  on  either  the  (lll)A  or  (100)  faces  of  InP 
substrate  crystals.  This  was  concluded  primarily  because  of  the  ease  in 
preparing  samples  suitable  for  X-ray  analysis.  In  this  year's  effort, 
more  than  140  growth  runs  have  been  made  employing  more  than  150  substrates. 
Of  these,  over  45  rims  were  made  on  (lll)A,  14  on  (100),  and  the  remainder 
on  (lll)B  oriented  substrates.  All  the  effort  was  directed  primarily 
toward  the  composition  Gag  27In0  73^S0  60P0  40  althou8h  other  compositions 
have  also  been  grown.  On  (lll)B  substrates  this  composition  was  realized 
last  year  and  can  now  be  grown  reproducibly  and  routinely.  On  the  (lll)A 
face,  however,  the  phase  diagram  is  unknown  in  this  composition  range. 

Hence  the  epilayer  compositions  obtained  encompassed  a  relatively  wide 
range  due  to  the  exploratory  nature  of  the  investigation. 

2.2  Phase  Diagram  Data 

During  liquid  phase  epitaxy,  growth  parameters  such  as  distribution 

coefficients,  k^,  for  different  atomic  species  (k^  =  C^/C^,  w^ere  C®  and 
l 

are  concentrations  in  the  solid  and  liquid  phases,  respectively;  i  * 

Ga,  In,  As,  or  P)  and  layer  morphology  are  dependent  upon  substrate  orien¬ 
tation.  Therefore,  to  obtain  the  same  epitaxial  layer  composition  on 
different  substrates,  the  melt  composition  has  to  be  varied.  Thus  liquidus 
isotherms  of  the  phase  diagrams  have  to  be  determined  for  new  melt 
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composition  ranges.  The  experimental  procedure  used  for  this  is  the  source 
dissolution  method  described  in  Reference  7.  Briefly,  a  melt  with  known  Ga 
atom  fraction,  x^a  (or  the  In/Ga  ratio),  is  prepared,  and  a  known  amount  of 
As  (in  the  form  of  InAs)  is  added.  Then  this  melt  is  saturated  with  an  InP 
source  wafer  to  establish  thermodynamic  equilibrium  between  InP  and  the 
melt.  It  has  been  found  that  a  period  of  two  hours  at  650° C  is  adequate 
for  this  process  for  small  (3  to  4  g)  melts.  From  weight  loss  measurements 
on  the  InP  source  wafer,  equilibrium  melt  composition  can  then  be  deter¬ 


mined.  Numerous  such  experiments  were  conducted,  and  the  results  are  shown 


in  Figure  1.  Here  it  is  seen  that  for  any  given  gallium  concentration, 

£ 

decreases  with  increasing  x^,  as  is  to  be  expected.  Also  it  is  to  be 


i 

x 

P 


noted  that  only  a  small  amount  of  phosphorous  is  present  in  the  melt. 


reflecting  its  high  distribution  coefficient  (k  ).  Further,  the  small 

P 

amount  of  InP  dissolved  in  the  melt  presents  problems  in  weight  loss  mea¬ 
surement.  In  Figure  2,  k.  is  plotted  as  a  function  of  y  for  Ga,  In  P.  As 

i  1— x  x  1-y  y 

epitaxial  layers  grown  on  (111) A  face.  Data  for  growth  on  (100)  and  (lll)B 


orientations  are  also  included  in  this  figure.  The  data  for  (100)  face 
were  obtained  from  published  work  [4,5]  and  our  own  work. 


2.3  Doping 

Three  donor-type  impurities — Sn,  Te,  and  Se — were  used  to  obtain 
.  17  -3 

carrier  densities  of  about  10  cm  or  above  in  the  Ga.  In.  P,  As  alloys. 

1-x  x  1-y  y 

All  three  were  found  to  be  well  behaved,  shallow-level  donors  with  a  doping 
behavior  similar  to  that  experienced  with  InP  but  with  slightly  different 
distribution  coefficients.  Our  data  are  presented  in  Table  1  for  alloy 
layers  grown  on  (lll)B  oriented  substrates  at  650°C.  Also  shown  are  data 
for  InP  [8,9]  for  purposes  of  comparison. 
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Figure  1.  Liquidus  curves  at  650 °C  showing  dependence  of  phosphorous 
solubility  on  concentration  of  arsenic  in  GalnAsP  melts. 
Running  parameter  is  gallium  atom  fraction. 


DISTRIBUTION  COEFFICIENT 


Table  1.  Distribution  Coefficients  of  Donor  Atoms  at  650°C 

Donor 

InP 

Ga0 . 27 In0 . 7  3As0 . 60P0 . 40 

Ge 

0.011 

Sn 

0.0022 

0.002 

Te 

0.4 

0.35 

Se 

4.2 

~1.0 

Si 

40.0 

2.4  Morphology 

Although  it  is  easy  to  chemically  process  (100)  oriented  substrates 
for  preparing  X-ray  analysis  samples,  it  has  other  limitations.  The  main 
problem  is  in  the  growth  of  thick  (25  to  30  pm)  layers.  When  layers  of 
such  high  thickness  are  precipitated  from  a  melt,  the  layer  composition 
changes  over  its  thickness,  resulting  in  a  changing  lattice  parameter. 

This  leads  to  a  severely  strained,  and  in  extreme  cases  polycrystalline, 
layer,  since  (100)  orientation  in  zinc-blende  structures  cannot  accommo¬ 
date  significant  lattice  mismatch.  This  problem  deserves  special 
attention  in  any  future  work. 

Generally,  the  crystalline  quality  and  morphology  of  thinner  (1  to 
5  pm)  layers  has  been  excellent  and  has  been  improved  substantially  over 
that  from  last  year.  This  is  evidenced  by  the  sharp  line  widths  obtained 
from  both  X-ray  diffraction  and  photoluminescence  peak  widths.  In  particu 
lar,  lattice  mismatch  between  epilayer  and  substrate  has  been  reduced  to 
below  the  limit  of  detection,  which  in  our  case  is  about  ±0.03  percent. 

The  control  of  layer  thickness  has  been  improved  and  is  better  than 
±20  percent  for  layers  1  to  5  pm  thick.  For  thinner,  submicron  layers, 
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the  growth  rate  variation  was  much  larger.  Further  work  is  needed  to 
achieve  precise  control. 


19. 
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3.0  MATERIAL  PROPERTIES 


3.1  Composition  and  Structure 

Composition  of  the  quaternary  epitaxial  layers  was  determined  by 
electron  probe  microanalysis  (EPMA).  In  addition,  single-crystal  X-ray 
diffraction  (XRD)  and  photoluminescence  (PL)  measurements  were  used  to 


measure  the  lattice  parameter  (a  )  and  energy  bandgap  (E  ),  respectively. 

o  g 


As  in  the  past,  aQ  and  E^  have  also  been  used  to  compute  the  composition 


of  epitaxial  layers  by  using  published  parametric  equations  (10,  11,  12], 
These  measurements  are  also  used  as  a  measure  of  the  quality  of  epitaxial 
layers:  Narrow  XRD  and  PL  peaks  indicate  high  crystalline  quality  and 
homogeneity.  Figure  3  illustrates  some  XRD  profiles  for  samples  with  zero 
and  nonzero  lattice  mismatch.  In  this  figure  the  intensity  of  the  dif¬ 
fracted  X-ray  beam  is  plotted  (in  arbitrary  units)  against  the  diffraction 
angle,  26.  A  Cu  target  tube  was  used  for  this  measurement,  which  yielded 


a  doublet  pair,  Ko^-Kot^,  with  wavelengths  of  1.5405  and  1.5443  A,  respec¬ 
tively.  Bare  InP  substrate  crystals  of  good  crystalline  quality  can 
resolve  this  doublet  extremely  well  with  a  2  0  separation  of  0.6  degree 
and  a  FWHM  of  approximately  9'.  With  GalnAsP  quaternary  layers  on  InP 
substrates,  two  sets  of  diffraction  peaks  (Ka^,  Ko^)  are  distinguishable, 
arising  from  the  epitaxial  layer  and  the  substrate.  If  the  epitaxial 
layer  is  lattice  matched,  the  peaks  merge  together  and  only  one  pair  of 
peaks  is  visible.  In  Figure  3,  Sample  Q-2-016  shows  this  situation.  Here 
the  Ka^  peak  has  a  FWHM  of  approximately  12',  slightly  greater  than  the 
9'  obtained  for  the  underlying  InP  substrate.  Diffraction  profiles  for 
Samples  Q-2-015  and  Q-2-027  show  small  shoulders  to  the  left  and  right, 
respectively.  For  Q-2-016,  since  the  peaks  are  merged  completely,  it  can 
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only  be  deduced  that  the  mismatch  is  less  than  ±0.03  percent,  the  detection 
limit.  Exactly  what  is  represented  by  this  peak  broadening  is  not  clear 
yet,  although  minor  gradients  in  lattice  parameter  or  a  defect-strained 
lattice  could  certainly  be  responsible. 

Table  2  lists  PL  peak  widths  for  some  samples.  It  is  seen  that  these 
peaks  have  a  FWHM  ranging  from  about  50  to  80  meV  at  room  temperature, 
which  compares  favorably  with  published  data  [13]. 

As  reported  previously  [7],  X-ray  diffraction  measurements  had  revealed 
the  presence  of  weak,  but  detectable,  lines  which  signaled  the  presence  of 
modulations  in  the  lattice  constant  of  the  quaternary  alloy  sample.  These 
were  interpreted  as  possibly  due  to  some  tendency  toward  long-range  order¬ 
ing  in  the  crystal.  An  intensive  investigation  was  carried  out  in  an  effort 
to  see  whether  or  not  a  definite  determination  could  be  made. 

Very  high  sensitivity  step  scans  (using  quartz  crystal  monochomatized 

MoK  radiation)  were  taken  through  the  background  between  the  Bragg  peaks 
a 

along  the  reciprocal  lattice  lines  OOO-OhO,  OOO-hhO,  and  000-hhh.  Stability 
of  the  radiation  source  was  sufficient  to  ensure  reproducible  output  over  a 
time  period  of  many  weeks.  Care  was  taken  to  eliminate  fractional 

Table  2.  Room  Temperature  Photoluminescence  Peak  Widths 
_  .  Bandgap  FWHM  Mismatch 


(eV) _ (meV) _ (%) 


2-021 

0.99 

90 

* 

2-022 

0.99 

73 

-0.06 

2-024 

0.91 

82 

* 

2-031 

0.89 

64 

* 

2-035 

0.87 

54 

* 

2-040 

0.77 

105 

-0.03 

2-045 

0.86 

54 

0.24 

5-015 

0.81 

90 

0.26 

5-017 

0.86 

60 

0.24 

5-027 

0.74 

94 

0.30 

* 

Mismatch  <  ±0.03% 
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wavelengths  from  the  monochromator  by  properly  adjusting  the  pulse  height 
analyzer.  Operating  parameters  are  listed  In  Table  3,  and  diffraction 
patterns  are  shown  in  Figures  4  through  7.  Approximately  three  minutes 
were  spent  collecting  counts  for  each  angular  increment  shown. 

Generally,  diffraction  pattern  backgrounds  were  higher  in  all  direc¬ 
tions  of  the  recpirocal  lattice  in  the  quaternary  crystal  than  they  were 
in  that  for  the  InP  crystal  substrate.  However,  as  reported  previously, 
the  quaternary  alloy  yielded  a  number  of  well-defined  but  weak  Bragg 
peaks  where  the  intensities  were  three  or  more  orders  of  magnitude  below 
the  principal  Bragg  peak  intensity.  Weak  reflections  not  permitted  in  the 
zinc  blende  structure  were  observed  to  be  as  follows:  (110),  *5(111), 

*5(500)  ,  *s(550) ,  (330)  .  It  was  concluded  that  these  were  not  induced  by  frac¬ 
tional  wavelengths  since  not  only  was  the  possibility  of  having  fractional 
wavelengths  eliminated  by  use  of  a  carefully  adjusted  pulse  height  analyzer 
but  also  the  excitation  potential  was  selected  so  as  to  locate  X/2  at  the 
cut-off  wavelength  of  the  continuous  spectrum.  Furthermore,  the  absence 
of  extra  reflections  in  a  control  crystal  of  InP  verified  the  absence  of 
fractional  wavelengths.  Any  possible  double  reflection  resulting  from  com¬ 
bining  two  existing  reflections  was  excluded  by  the  extinction  rule  for 
zinc  blende  structure  where  h  +  k,  k  +  £.,  (£.  +  h)  =  2  n  for  (hkJt)  must  be 
satisfied.  Thus  it  was  concluded  that  the  extra  peaks  seen  were  due  to 
the  presence  of  some  structure  present  in  the  Ql-72  quaternary  crystal  in 
addition  to  the  normal  zinc  blende  structure.  Following  these  observations, 
a  variety  of  crystal  structures  was  systematically  examined  to  see  whether 
or  not  the  extra  lines  could  be  explained  by  any  other  structure. 

Speaking  in  general  terms,  long-range  ordering  results  in  a  derivative 
structure  from  the  basic  structure  with  a  set  of  equipoints  being  divided 
into  two  or  more  subgroups  (with  an  increase  of  unit  cell  size  if  necessary) . 
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Pulse  height  analyzer  setting  Upper  level:  0.8  keV 


200  fOO  600  800 
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Figure  7.  29-Step  scan  pattern  showing  extra  reflection  %  (111). 
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The  zinc  blende  structure,  which  is  the  basic  structure  of  the  derived 
superlattice,  belongs  to  the  cubic  system.  To  check  crystal  symmetry, 
equivalent  reflections  IlO,  101,  and  Oil  were  chosen  for  their  identical 
absorption  corrections.  Three  peaks  of  identical  shape  and  size  were 
observed  as  shown  in  Figures  4,  8,  and  9,  showing  that  (111)  is  a  threefold 
axis  and  leading  to  the  conclusion  that  the  structure  is  either  cubic  or 
rhombohedral . 

By  doubling  the  edge  length  of  the  unit  cell,  the. observed  superstruc¬ 
ture  reflections  could  be  indexed  properly  and  become  220,  111,  500,  550, 
and  660  of  the  larger  cell.  Thus  it  was  concluded  that  there  is  a  distinct 
possibility  of  long-range  ordering  on  the  anion  or  cation  sublattice,  or  on 
both.  Further,  it  can  be  concluded  that  the  cubic  unit  cell  of  such  a 
superlattice  contains  64  atoms  rather  than  eight  and  that  the  unit  cell  has 
a  lattice  constant  twice  that  of  the  host  lattice. 

Since  the  coordinates  for  all  the  atoms  in  the  unit  cell  are  known, 
the  usual  search  for  the  possible  space  groups  implied  by  "diffraction 
symbol"  through  determining  conditions  for  a  set  of  possible  reflections 
is  unnecessary  and  time  consuming.  A  survey  of  the  36  possible  cubic  space 
groups  capable  of  allowing  correct  placement  of  cations  and  anions  deter¬ 
mined  that  the  possible  space  groups  are  limited  to  the  following:  P23, 
P2l3,  F23,  123,  and  P^3n.  Coordinates  for  the  special  and  general  posi¬ 
tions  of  each  of  the  equipoints  are  listed  in  Table  4.  The  observed 
superstructure  reflections  111,  200,  500,  550,  and  660  satisfy  only  space 
group  P23  and  none  of  the  others. 

Further  discussion  of  structural  possibilities  is  provided  in 
Appendix  A. 
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Figure  8.  20-Step  scan  pattern  showing  extra  reflection  (loi) . 


Possible  Space  Groups  for  Cubic  Structure  Model 


•  X  •— *|  00 


l' 

ill 


•  X  cn|co 


X  en|oo 


0) 


•-Mcm 

•—•I  •4' 

>. 

1  Sk 

»  >, 

>4 

o 

a 

X 

•  X 

X 

l  X 

* 

N 

N 

1  N 

1  N 

O  i-4|CM 

•* 

N 

1  N 

1  N 

N 

O  ^|<M 

n 

>» 

1  >• 

>> 

1  * 

«  a> 


io|  co  m|o  x  n  js  <o|co 

g  g  > 

XX  XXX 


N 

,  N  IN 

1  1  _ 

CO  £>  ^H|CM  T*  .O 


1-4 

+ 

•5k 

1-4 1  CM 

*-l 

i-4)cm 

5k 

>» 

X 

»  5k 

i 

+ 

+ 

9k 

X 

r-l|<M  i-4|cM 

i-4|CM 

i 

* 

• 

IX 

•k 

X 

4fc 

i  X 

o 

1-4|CM 

4k 

X 

+ 

i  X 

•k 

g*"N 

H|ki 

«k 

X 

g 

N 

.  N 
+ 

«-4l<M 

«k 

N 

1 

0 

N 

■ 

— 1| 

1-4 | CM 

t  N 

r-t|04 

1 

-4|cm 

% 

1  N 

N 

N 

•* 

9k 

X 

X 

•“•I'* 

«k 

1 

.  + 

O 

+ 

*-4|cN 

i 

•— l|CM 

g 

fs 

Jk 

1 

*-*|CM  *H|c>4 

•k  «k 

1 

>> 

«k 

t-4|CM 

>4 

l>» 

X 

1 

/■—S 

•k 

9k 

+ 

«k 

O  X 

T 

X 

1-4 1 04 

X 

col* 

a  x 

•k 

u 

II 

+ 

«k 

1 

g 

X  X 

l  X 

^x 

X 

4—4 1  CM 

i  X 

1-4 1  CM 

X 

N— -> 

/^k 

n 

9 

■w 

•o 

w 

e 

<■ 

<* 

04 

1*4 

04 

H 

! 


CO  H 
0\  04 


t-<lCv|  *"*l  CM  rH|CM 

*  M  •  « 

X  X  X  X 

+  *  +  «  ^Too 

W|<vj  ,-l|c\i 


X  X  X  X 
+  +  •  • 

—•I  CM  r-»|  e>l  «— l|CM  »H|eM 


XX 

J  !  «*M°0 


•  X  «x 
l  X  X 
X  I  X 


O 

•-•I  C4  o 

r-*I  <N  O 


X  •  X  I  X  X 

N  N  I  N  I  N 
IJS  >\  •>> 


h|n  —M  cv 


N  ►> 

+  » 

•*  •» 

X  N 

+  I 

»-i|cm  i-MC'J 

*  ft 

>»  X 

+  + 

•— *  |  OvI 

JS  X 

+  + 

*-<|«S  *-l|cM 


N  >»  X 

I  +  + 

— *|  CVi  r-4|cM  >-l|  CM 

ft  ft  ft 

X  N 

+  I  I 

<-H|cM  r— 1|  CV4  r-(|  CS 


*— 1|  CM  r-«|CM  f-l|  CS1 


*— 1|  CM  r-l|CM  —l|0< 


**(<M  ©  «-<|CM 

«-»|?N  •— *1  <N  O 

a 

X 

N 

X 

1  X 

l  X 
N 

l  X 
l  N 

N 

+ 

>> 

+ 

X 

1 

N 

+ 

>* 

+ 

•  X 

i 

>-l|  CN 

I-I|CM 

* 

«k 

•ft 

•ft 

*■ 

•ft 

O  H|N  O 

N 

1  N 

1  N 

N 

X 

N 

X 

N 

>> 

O  H|N  O 

i 

►> 

■  >» 

+ 

+ 

+ 

1 

1 

1 

©  © 

X 

X 

•  X 

X 

IX 

*-l|CM 

i—l  t  CM 

i-H|CS 

1-<|CS 

*H|CM 

«  JS 


CM  SO 


«o  c 

r-t  n 
M  Kf 


mmnrnmmm 


At  this  point  the  existence  of  long-range  order  has  not  been  unequivo¬ 
cally  established.  In  fact,  as  discussed  in  Section  5,  there  is  evidence 


that  at  least  some  of  the  quaternary  alloy  samples  prepared  in  this 
research,  if  not  all,  have  considerable  defects,  the  nature  of  which  has 
yet  to  be  determined. 

3.2  Electrical  Properties 

Hall  effect  measurements  were  made  using  the  standard  Van  der  Pauw 

method  to  determine  Hall  mobility  and  net  carrier  concentration.  For 

unintentionally  doped  layers,  the  net  electron  concentration  due  to  back- 

16  “3 

ground  impurities  is  found  to  be  in  the  low  10  cm  range.  A  summary  of 
these  measurements  is  given  in  Table  5.  For  higher  doping  density,  three 
donor  impurities — Sn,  Te,  and  Se — have  been  used  successfully.  All  three 
are  well  behaved,  shallow  donors  in  this  quaternary  with  a  doping  behavior 
similar  to  that  in  InP  but  with  slightly  different  distribution  coefficients, 
as  discussed  previously  in  Section  2.3. 


Table  5.  Carrier  Concentrations  and  Electron  Mobility  Values 
at  300  K  and  77  K  for  Quaternary  Alloy  Samples 


-17  -3 

n(x  10  cm  ) 

,  2 
PQ(cro 

/V-s) 

pie 

300  K 

77  K 

300  K 

77  K 

22 

2.75 

—  ^  — 

2248 

... 

24 

0.114 

0.132 

2282 

4990 

25 

0.103 

— 

1952 

— 

32 

0.762 

0.734 

3174 

4295 

35 

3.43 

3.16 

3074 

4086 

10 

0.483 

— 

1368 

— 

26 

1.64 

— 

2071 

» 
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Using  the  form 


Che  samples  for  which  300  K  and  77  K  data  are  reported  in  Table  5  show  a 
temperature  dependence  which  ranges  from  b  ■  -0.21  to  b  =  -0.58.  As 
expected,  this  is  a  definite  correlation  between  the  carrier  concentration 
and  the  reduction  in  mobility. 

Theoretical  calculations  of  electron  mobility  values  in  Ga,  In  P,  As 

1— x  x  1-y  y 

alloys  using  Monte  Carlo  techniques  have  been  made  for  those  alloy  composi¬ 
tions  which  are  lattice  matched  to  InP  [6].  The  theoretically  predicted 
variation  of  Hall  mobility  is  shown  in  Figure  10  compared  to  the  range  of 
values  obtained  experimentally  [14,15].  Efforts  to  achieve  a  better  fit  of 
theory  to  experiment  [16]  have  had  to  invoke  relatively  large  amounts  of 
compensation  by  acceptor  impurities  to  increase  the  contribution  of  ionized 
impurity  scattering  as  well  as  large  amounts  of  "alloy"  scattering.  The 

latter  has  been  treated  theoretically  by  several  approaches  [17]  which  give 
1/2 

a  T  dependence — namely,  so-called  "space-charge  scattering"  [18,19] — 

which  acts  to  decrease  the  mobility  below  that  predicted  for  III-V  compound 
alloys  when  the  mechanisms  of  polar  optical  mode,  ionized  impurity,  and 
alloy  scattering  are  invoked. 

Figures  11  and  12  show  data  points  derived  from  Hall  effect  and  conduc¬ 
tivity  measurements  by  Houston  and  his  co-workers  [20]  on  two 
Ga^_xInxPj_yASy  alloy  samples  grown  by  standard  liquid  phase  epitaxy  tech¬ 
niques  in  our  laboratory.  These  layers  were  grown  on  iron-doped  InP 
substrates  and  are  both  almost  exactly  lattice  matched  to  InP.  The  Hall 
mobility  for  both  samples  rises  to  a  maximum  at  about  100  K  in  both  cases 
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□  Greene  et  at.  (1979) 

6  Houston  t  at.  (1978) 
O  RTI/NCSU 
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Figure  10. 

Comparison  of  theoretical  and  experimental  values  for 
mobility  in  Ga.  In  P.  As  epitaxial  layers  lattice 
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and  then  decreases  to  a  saturating  value  at  temperatures  below  about  10  to 
20  K.  The  latter,  of  course,  is  explicable  in  terms  of  impurity  band  con¬ 
duction  [21].  The  shape  of  the  experimental  curves  in  Figures  11  and  12 
may  be  explained  by  the  combination  of  polar  optical  mode,  ionized 
impurity,  and  alloy/space  charge  scattering  alone  or  by  invoking  an  addi¬ 
tional  mechanism,  dislocation  scattering  [22].  A  reasonable  fit  to  the 
experimental  curves  is  obtained  by  using  a  simple  extension  of  the  models 
used  for  ternary  alloys  [23]  and  adapting  the  Dexter  and  Seitz  [22] 


formula. 


3.05  x  103  EV/i^Liv\2  *±(SL.\(V^_' 
1  \1'Vj  '  WA  cm2 


where  is  the  deformation  potential  (eV) ,  X  is  the  crystallographic  slip 

distance  (cm),  v  is  the  Poisson  ratio  for  the  material,  Nd  is  the  number 

2 

of  dislocations  per  cm  ,  and  T  is  the  temperature  (K) .  Nd  was  treated  as 

an  adjustable  parameter,  and  estimates  were  used  for  the  other  parameters 

—8 

as  follows:  E^  **  6.0  eV,  X  =  2.5  x  10  cm,  and  v  =  0.3. 

Weisberg's  formula  [19]  for  space  charge  mobility  was  used  in  the  form 


1  ,  „  --10  / m*\1/2  T1/2  /v  -  s 

^  \m0/  X  V  cm2  . 


where  the  parameter  X  (in  cm)  is  the  mean  free  path  for  scattering  by  large 
cross  section  objects  (opaque  scattering  centers)  whose  cross  section  is 
taken  in  this  treatment  to  be  independent  of  temperature.  This  mean  free 
path  is  equal  to  the  product  of  N  ,  the  number  of  opaque  scattering  centers 
per  unit  volume,  and  A^,  the  cross-sectional  area  of  the  scattering  center. 
More  properly. 
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“s 

E  A 

k=l 


Sk 


where 


1  S 

As  =  n7  ,E,  ASk 


S  k=l 


(4) 


(5) 


and  A  represents  the  cross-sectional  area  for  electron  scattering  for 

u  tv 

the  k^  object  in  the  unit  volume. 

Fermi-Dirac  integrals  were  used  to  calculate  energy  averages  for  the 
mobility  and  scattering  factors  for  each  mechanism.  For  Curves  A  and  B 
for  each  sample  the  effect  of  degeneracy  on  the  screening  radius  for  ionized 
impurity  scattering  has  been  taken  into  account  [24] . 

Combination  of  reciprocal  mobilities  for  the  various  scattering  mecha¬ 
nisms  is  used  to  calculate  a  conductivity  mobility.  This  is  then  corrected 
by  an  average  scattering  factor  calculated  from 

£  V^V 

-  (6) 

<1/v 

k  K 

to  obtain  the  Hall  mobility. 


(7) 


As  can  be  seen  from  Figures  11  and  12,  good  fits  are  obtained  in  the  tem¬ 
perature  range  50  to  300  K  for  reasonable  values  of  the  parameters. 

—6  0 
For  both  samples  A  ~  6.0  x  10  cm.  If  an  average  diameter  of  50  A 

17  -3 

is  assumed,  Ng  ~  8  x  10  cm  .  In  this  case  approximately  one  half  of 
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Che  volume  would  be  occupied  by  Che  hypoCheCical  opaque  scaCtering  centers. 
Since  Che  fracCion  of  volume  occupied  by  the  (assumed  spherical)  opaque 
scattering  centers  is  given  by 


fs  -  f  dsrt  ,  (8) 

o 

reduction  of  the  diameter,  dg,  to  10  A  would  reduce  the  volume  fraction, 
f  ,  to  about  9  percent. 

If  the  opaque  scattering  centers  occupied  an  appreciable  fraction  of 

the  volume  of  the  alloy  and  were  of  different  crystal  structure,  one  would 

expect  detection  by  X-ray  diffraction.  However,  our  X-ray  diffraction 

measurements  of  similar  Ga  In  P  As  alloy  samples  prepared  in  our 

x  x  x  x  y  y 

laboratories  have  not  revealed  the  presence  of  any  lines  that  could  be 
interpreted  as  the  existence  of  a  separate  phase  as  inclusions  (e.g., 
indium).  As  discussed  previously  in  Section  3.1,  extra  lines  do  appear  [25], 
but  these  exhibit  the  same  symmetry  as  the  host  lattice  with  an  apparent 
doubling  of  the  lattice  constant.  Therefore,  either  the  volume  fraction  is 
small  enough  to  reduce  diffraction  lines  down  to  the  noise  level,  or  the 
opaque  scattering  centers  have  essentially  the  same  symmetry  properties  as 
the  host  lattice. 

Since  there  is  some  controversy  [17]  about  the  proper  choice  for  the 
scattering  potential  used  In  the  calculation  of  alloy  scattering,  use  of 
differences  in  energy  bandgap  values  rather  than  either  electron  affinity 
differences  or  electronegativity  differences  was  examined.  This  would 
increase  the  alloy  scattering  rate  by  a  factor  of  about  three  whereas  an 
increase  by  a  factor  of  about  24  is  required  to  explain  the  difference 
necessary  to  fit  the  experimental  data. 
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There  is  evidence  that  the  effects  producing  the  shape  of  the  mobility 

curves  for  the  Ga,  In  P.  As  are  caused  by  some  mechanism  other  than 
1-x  x  1-y  y 

alloy  scattering  per  se.  Figure  2  in  Reference  26  displays  a  plot  of 

electron  mobility  versus  temperature  for  vapor  phase  epitaxial  InP  samples. 

One  of  these  samples,  CV66,  with  a  nominal  electron  concentration  of 
17  -3 

3.5  x  10  cm  at  300  K,  has  a  curve  shape  similar  to  those  for  the  alloy 
sample  shown  in  Figures  11  and  12  of  this  report.  Similarly,  Figure  3  in 
Reference  27  displays  a  plot  of  mobility  versus  temperature  for  samples  of 
GaAs  obtained  commercially  £a  1970  (presumably  either  Czochralski  or 
Bridgman  process  growth) .  Two  of  these  curves  have  shapes  similar  to 
those  for  the  alloy  samples  in  Figures  11  and  12. 

In  summary,  for  many  materials  with  moderate  values  of  electron  con¬ 
centration  (~10'*'6  to  lO1^  cm  ,  there  appears  to  be  a  mechanism  present 

that  contributes  a  component  of  scattering  in  the  range  100  to  300  K  at  a 

1/2 

rate  proportional  to  T  .In  alloys,  this  process  is  indistinguishable 

from  so-called  alloy  scattering  except  for  those  cases  in  which  the  magni- 

1/2 

tude  of  the  T  component  of  scattering  considerably  exceeds  that 
calculated  from  the  simple  theory  for  alloy  scattering. 

In  the  range  of  about  50  to  100  K  the  mobility  in  the  samples  shown 
in  Figures  11  and  12  appears  to  be  controlled  either  by  ionized  impurity 
scattering  or  by  a  combination  of  this  with  scattering  from  dislocations. 
The  dislocation  density  calculated  on  the  basis  of  Eq.  (2)  and  the  assump¬ 
tions  made  concerning  the  parameter  values  used  in  this  equation  are  in  the 
10  -2 

range  of  2  x  10  cm  .  Both  samples  were  almost  exactly  lattice  matched 

to  the  InP  substrates,  and  the  substrates  themselves  had  etch  pit  densities 
_2 

of  about  6,000  cm  .  Thus  the  existence  of  such  a  large  dislocation 
density  cannot  be  attributed  to  either  of  these  factors.  However,  if  the 
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opaque  scattering  centers  invoked  previously  are  due  to  the  presence  of  a 

phase  with  a  lattice  constant  sufficiently  different  from  that  of  the  host 

lattice,  such  high  densities  are  plausible.  Using  a  simple  calculation  [28], 

10  2 

the  lattice  mismatch  required  to  generate  2  *  10  dislocations  per  cm  is 
estimated  at  about  0.8  percent. 

Based  on  the  evidence  presented,  it  appears  that  the  electron  mobility 
in  Ga^  ^In^P^  As  alloys  grown  by  liquid  phase  epitaxy  may  be  severely 
reduced  by  the  presence  of  one  or  more  phases  that  act  both  as  opaque  scat¬ 
tering  centers  and  may  in  addition  act  as  generators  for  dislocation 
networks.  However,  this  phenomenon  is  also  observed  in  InP  grown  by  vapor 
phase  epitaxy  [26]  and  by  bulk  growth  GaAs  [27].  Since  both  of  these 
binaries  are  also  obtained  with  high  quality  and  excellent  electronic 
properties,  it  is  concluded  that  there  is  considerable  hope  for  improving 
the  quality  of  Ga^^In^^^As^  alloys. 


4 . 0  DEVICE  TECHNOLOGY 


4.1  Contacts 

Ohmic  contacts  to  quaternary  layers  have  been  made  using  evaporated 
In  or  Au-Sn  contacts  with  subsequent  alloying  at  400°C  for  three  minutes 
in  a  hydrogen  ambient. 

For  forming  Schottky  barriers,  several  metallization  processes  have 

0 

been  used — e.g.,  Au,  Al,  Cr,  Ag,  Mg.  Thin  layers  (~1,000  A)  of  these 
metals  were  deposited  on  precleaned  samples  through  a  metal  shadow  mask 
with  circular  openings.  Precleaning  was  done  with  a  thorough  rinse  in 
organic  solvents  followed  by  a  brief  15-second  etch  with  0.1-percent 
bromine -methanol  solution  with  a  final  methanol  rinse.  After  many  trials, 
we  have  finally  been  able  to  get  rectification  with  Au.  Figure  13  shows  a 
C-V  plot  for  a  203- pm  diameter  metal  dot.  Here  the  barrier  height  is  seen 

to  be  ~0.51  V  for  a  quaternary  material  bandgap  of  0.96  eV  ( d>  /E  =  0.54). 

B  g 

In  Figure  14  the  forward  I-V  characteristic  of  this  device  is  shown.  From 
a  straight  line  fit  in  the  low  current  injection  range,  it  is  determined 
that  <j>  =  0.51  V  and  the  diode  factor  n  *  2.6. 

D 

Attempts  have  been  made  to  make  Schottky  barriers  on  other  quaternary 
sanples  with  Au  metal,  but  all  these  have  proven  to  be  unsuccessful.  The 
reasons  for  this  behavior  are  not  yet  fully  understood.  From  all  the 
experiments  that  have  been  performed  in  this  period,  no  systematic  trend 
for  the  occurrence  of  a  rectifying  barrier  is  evident.  It  is  clear  that 
more  work  is  required  in  at  least  two  areas:  (1)  in  preparation  of  clean 
surfaces  and  (2)  in  structural  evaluation  of  quaternary  material  for  defect 


structures. 


Figure  13.  C-V  curve  for  Au/n-GalnAsP  Schottky  barrier  diode. 


CURRENT  (A) 
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4.2  MESFET  Design 

Preliminary  designs  undertaken  for  MESFET  devices  optimized  for  X-band 
(8  to  12  GHz)  and  K  -band  (26  to  40  GHz)  frequency  ranges.  Details  are 

O 

given  in  Appendix  B.  A  summary  of  the  calculated  device  DC  parameters  and 
figures  of  merit  are  given  in  Tables  6  and  7,  respectively.  It  should  be 
pointed  out  that  these  results  were  based  on  theoretical  values  for  low 
field  mobility  and  peak  velocity  [Ref.  6],  As  reported  previously,  the 
values  of  low  field  mobility  obtained  to  date  have  been  about  half  that  pre¬ 
dicted  theoretically.  This  and  the  problems  cited  in  obtaining  good  Schottky 
barriers  have  prevented  device  fabrication  to  achieve  the  figure  of  merit 
values  in  Table  7. 


Table  6.  Device  DC  Parameters 


X-Band  Device 

K  -Band  Device 
a 

h"  <■» 

427 

136 

xds  (mA) 

64 

20.5 

Vds  <v> 

4 

2 

VDinch-off  (V) 

4.28 

1.64 

Table  7. 

Calculated  Figures  of 

Merit 

X-Band  Device 

K  -Band  Device 
a 

Cutoff  frequency,  f^ 

23.9  GHz 

95.7  GHz 

Max.  frequency  of 
oscillation,  f 

max 

51.2  GHz 

165  GHz 

Unilateral  gain,  U 

12.6  db  at  12  GHz 

12.3  db  at  40  GHz 
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I.  CUBIC  STRUCTURE  MODELS  AND  STRUCTURE  FACTOR  CALCULATIONS 


In  order  Co  determine  Che  exact  crystal  structure  of  the  super- 
lattice  formed,  it  was  necessary  to  set  up  various  trial  structure 
models,  and  then  calculate  the  structure  factors  for  each  of  these 
models.  Actual  observed  structure  factors  were  then  compared  with 
those  of  the  models. 

By  combining  the  equipoints  in  different  ways,  and  leaving  out 
those  that  might  be  correlated  by  a  simple  translation,  452  arrange¬ 
ments  compatible  with  the  273  hypothetical  compositions  are  obtained. 
This  i3  summarized  below. 


Ratio 

Permutation 

gy£.e~ 

for  Arrangemeni 

1 

aaaa 

4 

27 

abed 

8  or  16 

11 

abcc 

4  or  8  or  16 

3 

abab 

8  or  16 

1x4+  <13x16+14x8)  +  (2x4+7x8+2x16) 

+  (2x8+1x16)  *  452  .  . 

Since  superlattices  found  in  binary  system  are  mostly  in  simple 
ratio,  132  structure  models  (having  composition  ratios  up  to  3:5:4:4) 
were  set  up,  and  structure  factors  were  calculated  for  each  model. 

The  computer  program  XFLS4  (May  1977  version)  by  Levy  and  Busing 
was  used  for  structure  factor  calculation.  This  program  is  based  dn  an 


idea  of  calculating  through  symmetry  operations  and  thus  it  may  be 
applied  to  all  space  groups. 

Structure  factor  for  the  ideal  composition  ®a3^ni3^2^s14  ^or 
Gart  -_Inn  Q,Prt  ,-Asrt  Q_)  is  also  calculated,  which  is  the  nearest  type 
to  the  experimental  composition.  The  observed  intensity  ratio  of  super¬ 
structure  to  normal  peaks  tended  to  be  lower.  Due  to  the  associated 
thermal  and  compositional  disorder,  it  is  unknown  whether  it  might  be 
fully  ordered  in  this  ratio,  or  whether  it  might  be  stable  in  some 
simpler  ratio.  Among  the  132  structure  models,  some  of  the  calculated 
structure  factors  for  certain  reflections  are  zero,  although  they  have 
been  observed  in  experiment.  Where  this  is  the  case,  these  models  may 
readily  be  excluded.  Despite  this,  the  possible  number  of  arrangements 
still  amounts  to  hundreds . 

Based  on  these  primary  analyses,  the  diversified  complexity  of 
the  ordering  in  this  quaternary  crystal  is  conceivable.  It  is  thus 
proposed  that  further  study  first  be  limited  to  ternary  systems.  This 
would  allow  ordering  to  take  place  either  in  cation  or  in  anion  sub- 
lattice  sites  only,  and  the  ordering  then  degenerates  into  a  pseudo¬ 
binary  case.  In  order  to  test  the  validity  of  this  conclusion,  two 
ternary  compositions  were  step  scanned,  GaQ  q^°q  99Asq  6gPQ  34 
and  Ga^  ^In^  53^®  (Q2-012)  .  (It  should  be  pointed  out  to  the  reader 
that  Ql-75  has  only  0.01  percent  Ga  and  may  therefore  be  considered  to 
be  truly  ternary.) 
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II.  OTHER  STRUCTURE  MODELS 


Since  both  GaAs  and  InP  possess  zinc  blende  structure,  it 
is  expected  that  this  quaternary  compound  has  an  adamantine  (diamond¬ 
like)  structure  (necessary  conditions  as  electron  rules  summarized  by 
Parthe  [29]  are  satisfied).  The  adamantine  structures  form  a  subgroup 
of  the  tetrahedral  structures  and  show  geometrical  features  of  wiirtzite, 
zinc  blende,  and  related  structures. 

Parallel  to  the  analysis  based  on  zinc  blende  structure,  vari¬ 
ous  wiirtzite  derivative  structures,  known  polytypes  of  ZnS  ,  CuPt 
rhombohedral  structure,  and  chalcopyrite  structure,  are  tried  out  to 
coincide  with  the  observed  superlattice  reflections.  This  is  summar¬ 
ized  in  Table  8.  The  first  three  columns  represent  the  experimental 
reflections.  In  the  fourth  column,  the  reflections  that  coincide  with 
those  obtained  by  doubling  the  edge  length  of  the  original  cell  are 
listed.  One  can  find  the  corresponding  reflections  by  going  horizon¬ 
tally  in  the  table.  When  a  hexagonal  poly type  is  presumed,  this  can  be 
verified  by  a  loss  of  cubic  symmetry  in  the  experimental  crystal.  How¬ 
ever,  the  possibility  must  be  allowed  that  hexagonal  "shifts"  occur  on 
all  four  cubic  till}  planes,  which  might  restore  an  apparent  cubic  sym¬ 
metry.  The  same  rationale  can  be  used  to  invoke  the  tetragonal  chalco¬ 
pyrite  type  to  explain  extra  reflections.  If  the  tetragonal  c-axis 
alternates  between  all  three  cube  axes,  the  apparent  symmetry  of  this 
structure  would  be  cubic. 
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A  ternary  case  leading  to  a  primitive  cubic  cell  of  the  approxi¬ 
mate  size  of  GaAs  would  be  obtained  by  having  the  two  like  ion  species 
ordered  in  one  sublattice  as  in  Cu^Au  .  From  the  known  theory  of  that 
structure,  it  immediately  follows  that  all  primitive  cubic  reflections 
listed  in  the  column  below  3C  in  Table  8  should  be  observed.  This  is 
not  compatible  with  the  experimental  observations,  so  that  arrangement 
can  be  excluded  from  further  discussion. 

The  rhombohedral  CuPt  type  provides  an  ordering  scheme  that 
could  explain  the  doubling  of  the  cell  edge  of  GaAs  .  The  reflections 
due  to  this  atomic  arrangement  are  listed  under  CuPt  II  in  Table  8. 

It  cannot  explain  the  observed  reflections  in  the  experimental  crystals. 

The  orthorhombic  BeSiN^  structure  type  is  derived  from  the 
hexagonal  wiirtzite  type  by  making  a  *  2a^ex  ,  b  -  /  3  a^ex  ,  c  *  chex 
(referring  to  Parthe  [29]).  According  to  Table  8,  it  cannot  account 
for  all  of  the  observed  reflections  in  the  three  experimental  crystals. 

The  well-known  tetragonal  chalcopyrite  structure  is  obtained 

by  doubling  one  cell  edge  of  the  cubic  zinc  blende  type  so  that 

a  ■  a  .  .  and  c  =  2  a  .  .  .  It  could  explain  some  of  the  experi- 

cubic  cubic 

mentally  observed  reflections,  but  does  not  account  for  the  111  and  550 


reflections  of  the  larger  cubic  cell  with  a'  =  2  a 


GaAs 
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Bragg  Reflections  of  Conventional  GaAs  Derivative  Structure 
Coinciding  with  Primitive  Cubic  Reflections  for  a  =  11.30 


*Standard  cubic  reflections  of  the  GaAs  type  are  in  square  brackets. 
^Additional  reflections  for  the  primitive  structure  are  in  parentheses 
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III.  DISCUSSION  AND  RECOMMENDATION  FOR  FUTURE  RESEARCH 

The  present  approach  is  purely  geometrical  and  is  based  on 
crystal  symmetry  arguments.  No  detailed  energetic  consideration  is 
given,  since  it  is  as  difficult  to  theoretically  determine  which 
structure  is  most  stable  in  any  particular  case,  as  it  is  to  decide 
whether  a  given  metal  will  crystallize  in  a  FCC,  BCC,  or  HCP  crystal. 

In  both  cases,  the  criterion  is  the  maximization  of  the  cohesive 
energy. 

After  all  these  surveys,  the  superstructure  reflections  found 
in  two  ternary  and  one  quaternary  compositions,  are  not  explained  by 
any  simple  atomic  arrangement.  Structures  that  readily  suggest  them¬ 
selves,  namely,  the  ZnS  polytypes  2H,  4H,  6H,  8H,  9R,  10H,  12R,  15R, 

21R,  the  tetragonal  chalcopyrite  structure,  the  orthorhombic  BeSiN^ 
structure,  ordering  of  the  anion  or  cation  sub lattice,  as  in  Cu-Au  or 
Cu-Pt  ,  all  do  not  account  for  the  observed  and  missing  extra  reflec- 
t ions .  • 

The  experimental  observation  of  additional  reflections  in  the 

crystal,  of  conpos  trior.,  ^0. 19lnO.SlP0.10As0.90  •  Ga0.01ln0.99P0.34As0.66’ 
and  Ga^  ^In^  53^  suggest  that  some  derivative  superstructure  of  the 

GaAs  type  is  being  formed.  The  reported  tendency  toward  clustering  in 
Ga^  5InQ  -^P  [30]  seems  to  contradict  this  conclusion,  because  ordering 
and  clustering  represent  opposite  thermodynamic  tendencies  that  are  not 
usually  compatible.  It  should  be  pointed  out  that  the  present  study 
does  not  include  the  above  phosphide  and  chat  the  gal Hum- indium-arsenide 
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need  not  behave  as  the  phosphide.  On  the  other  hand,  it  is  also  pos¬ 
sible  that  the  alloys  studied  in  the  present  Investigation  segregate 
into  two  phases,  both  of  which  have  "tetrahedral"  structures.  The 
structure  of  the  new  phase  will  then  have  to  be  characterized  as  a 
derivative  of  GaAs  with  details  as  yet  undetermined. 

For  further  structure  analysis,  the  following  is  recommended: 

1.  Increasing  the  intensity  of  superstructure  reflections  is  essen¬ 
tial..  This  may  be  attained  by: 

(1)  Annealing  of  the  crystal  to  eliminate  the  thermal  disorder  factor. 

(2)  Using  Cuka  radiation. 

Shorter  wavelength  (Moka  radiation)  was  used  for  initial  study 
to  increase  the  number  of  accessible  points  in  reciprocal  lattice. 

Cuka  radiation  may  be  applied  for  further  study  for  its  imposed  effect 
of  modifying  atomic  scattering  factors.  This  is  estimated  as  the  fol¬ 
lowing: 

X  k 


Atomic 

No. 

(k  Absorb. 

Edge) 

Moka 

Vxk 

Cuka 

*/Xk 

Af 

In 

49 

0.444 

1.599 

3.474 

- 

Ga 

-  31 

1.195 

0.594 

1.290 

-2.2 

P 

15 

5.787 

0.123 

0.266 

- 

As 

33 

1.045 

0.679 

1.476 

-1.8 

X/Xk  < 

0.8,  A  f  correction 

is  negligible 

> 

•  X/Ak  > 

1.6,  A  f  correction 

is  about  constant 

and  Independent  of  variation 
in  Xk  . 
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The  contribution  of  Ga  and  As  scattering  to  those  reflections 
that  are  strongly  dependent  on  the  difference  f^a  -  can  be  en¬ 

hanced  by  using  Cuke  radiation  instead  of  Moka  radiation. 

loc(fr  -  f  )8  -  (31  -  33 )s  *  4 
«a  As 

lfcc(fGa  "  f As )2  “  1(31  -  2.2)-(33  -  1.8)]3  =  5.76  . 

2.  Rhombohedral  derivative  structure  model  should  be  applied.  It 
is  likely  the  case,  since  differential  thermal  expansion  along  rill] 
and  inclined  dislocations  [31,32]  resulting  from  LFE  growth  would 
account,  for  the  rhombohedral  structure,  according  to  Verna  and  Krishna 
[32].  Space  group  R3ra  (lowest  rhombohedral  symmetry  that  can  be  exhib¬ 
ited  by  a  periodic  stacking  of  close  packed  layers)  [32]  may  describe 
the  64-atom  unit  cell,  961  hypothetical  compositions ,  and  more  than 
one  thousand  arrangements  may  be  generated. 

Perhaps  the  best  and  quickest  clue  is  to  use  TEM  lattice  imaging 
for  primary  deduction  and  following  the  refining  by  x-ray  diffraction. 
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APPENDIX  B 

Proposed  Quaternary  Alloy  Semiconductor  MESFET  (QFET)  Designs 
for  Operation  in  X  and  Ka  Bands 


R.  J.  Trew 


I.  INTRODUCTION 


This  report  presents  proposed  device  designs  for  MESFETs  capable  of 
low-noise,  small-signal  amplification  of  RF  signals  in  the  X  (8-12  GHz) 
and  Ka  (26-40  GHz)  band  frequency  ranges.  The  devices  are  fabricated  from 
quaternary  alloy  Ga^  27In0  73P0  4As0  6  eP*tax*al  layers  grown  on  semi- 
insulating  InP  substrates.  Standard  one-micron  gate  length  technology  is 
to  be  utilized  for  the  fabrication  of  the  X-band  device.  Operation  in  Ka 
band  represents  a  significant  advance  in  the  current  state  of  the  art,  and 
improved  technology  will  be  required  in  order  to  fabricate  the  required 
submicron-sized  gate  length  devices.  The  Ka  band  structure  is  designed 
with  a  0.2  um  gate  length.  Standard  0.5-pm  gate  length  technology  could 
be  used  to  fabricate  the  Ka  band  MESFET  as  demonstrated  by  the  results  of 
Krumm  et  al.  [33].  Using  a  0.5-ym  GaAs  MESFET,  they  constructed  and 
operated  an  amplifier  for  Ka  band  operation  and  obtained  11  db  gain  at 
33  GHz  with  a  5.5-db  noise  figure.  A  shorter  gate  length  should  lead  to 
improved  performance. 

The  design  procedure  and  calculations  are  presented.  Design  compro¬ 
mises  were  made  on  the  basis  of  technological,  material,  and  device 
operation  criteria;  and  the  rationale  for  the  various  decisions  are  pre¬ 
sented.  Much  of  the  design  is  based  on  the  small-signal  model  of  Pucel, 
Haus,  and  Statz  [34].  This  model  has  been  used  by  various  workers  to 
accurately  predict  device  performance.  In  order  to  get  a  feeling  for  the 
performance  that  can  be  anticipated  from  these  devices,  certain  figures  of 
merit  are  calculated. 
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II.  DESIGN  PROCEDURE 


The  X  band  MESFET  is  shown  in  Figures  15-17  and  is  designed  with  a 
one-micron  long  gate.  The  material  and  design  parameters  are  listed  in 
Tables  9  and  10.  One-micron  gate  GaAs  MESFETs  are  extensively  used  in  the 
microwave  industry  for  a  wide  variety  of  applications  and  provide  excel¬ 
lent  performance.  With  the  same  technology  and  the  improved  parameters  of 
the  quaternary  alloy  semiconductor  material,  improved  device  characteris¬ 
tics  should  be  obtained.  The  frequency  response  and  noise  figure  of  the 
device  are  critically  dependent  upon  the  gate  length.  Although  the  one- 
micron  gate  is  near  optimum  for  X  band  operation,  shorter  gates  are  required 
for  Ka  band.  There  is  a  limit  to  performance  improvements  from  gate  length 
reduction  although  the  gate  length  limit  is  not  precisely  known.  Commer¬ 
cial  GaAs  MESFETs  are  available  with  0.5  ym  gate  lengths,  and  experimental 
devices  have  been  produced  with  0.25  ym  gates.  The  Ka  band  device  is 
designed  with  a  0.2  pm  gate,  and  state-of-the-art  fabrication  techniques 
will  be  required.  There  is  not  much  information  available  on  the  operation 
of  such  short  gate  length  devices,  and  the  RF  performance  may  be  signifi¬ 
cantly  affected  by  velocity  overshoot  considerations.  This  subject  requires 
further  research  the  Ka  band  MESFET  is  shown  in  Figures  18-20,  and  its 
design  parameters  are  listed  in  Table  10. 

The  width,  Z,  of  the  gate  affects  the  RF  output  power,  the  device  gain, 
and  almost  all  of  the  equivalent  circuit  elements.  It  has  a  significant 
effect  upon  the  gate  resistance,  Rg.  Although  the  gate  is  a  transmission 
line,  it  can  be  modeled  with  lumped  elements  as  long  as  the  gate  width  is 
short  compared  to  the  RF  wavelength.  In  order  to  avoid  distributed  effects, 
the  gate  length  should  be  no  more  than  approximately  0.04  A  in  the  long 
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Table  9.  Material  Parameters  for  Ga^  27In0  73P0  4As0  6 
Used  in  MESFET  Design 

CR  "  13‘2 

e  *  8. 854 (10) “14  F/cm 
o 

E  =  5 (10) 3  V/cm 
s 

v  =  2. 8(10) 7  cm/sec 
P 

v  (@  8  kV/cm)  *  2 (10) 2  cm/sec 

p  (x-band)  =  350  Q/D 
s 

p^(Ka-band)  =  366  fi/D 


Table  10.  Design  Parameters 


X-Band  Device  Ka-Band  Device 


Nd  (cm  3) 

1  x 

io17 

2.4 

x  1017 

-4 

Lg  (cm) 

1  x 

10  4 

0.2 

x  10 

-4 

-4 

Lsd  (cm) 

5  x 

10 

2.5 

x  10 

1.15 

-4 

Lgd  (cm) 

2.75  x 

10 

x  10 

-4 

a  (cm) 

0.25  x 

10 

0.1 

x  10 

-4 

100 

-4 

(cm) 

150  x 

10 

x  10 

-4 

-L 

(cm) 

300  x 

10  4 

100 

o 

t-i 

X 

N 

2 

2 
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Figure  19.  Ka-Band  MESFET — Gate  Detail. 
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dimension.  The  X  band  device  is  designed  with  a  300- um  gate  width,  and  the 
Ka  band  device  has  a  100-um  gate  width.  The  gates  of  these  devices,  there¬ 
fore,  can  be  modeled  as  lumped  elements.  The  gate  spacing  in  the  channel 
has  a  significant  effect  upon  the  operation  of  the  device  since  the  device 
capacitances  are  a  strong  function  of  geometry.  Generally  it  is  desirable 
to  minimize  the  gate-to-drain  feedback  capacitance,  and  the  gate  is  located 
asymmetrically  in  the  channel  with  the  gate-to-drain  spacing  larger  than 
the  gate-to-source  distance.  Also,  the  gate  channel  width  significantly 
affects  the  power  handling  capability  of  the  device  since  the  dc  voltage 
that  can  be  applied  is  dependent  upon  this  parameter.  The  limit  to  the  dc 
voltage  is  determined  by  the  drain- to-source  breakdown  voltage,  which  for 
GaAs  devices  can  be  expressed  approximately  as 

V  a  3  5  (i  \2/3 
fids  asd; 

For  high  power  devices,  therefore,  it  is  desirable  to  have  large  source-to- 
drain  spacings.  This,  of  course,  conflicts  with  high  frequency  response, 
since  the  device  is  transit-time  limited  at  high  frequencies.  The  X  band 
device  as  a  5-ym  spacing  with  the  gate  located  2.75  ym  from  the  drain.  The 
Ka  band  device  has  a  2.5-ym  gate  channel  with  the  gate  centered  in  the  chan¬ 
nel.  The  centered  channel  design  for  the  Ka  band  device  was  selected  to 
simplify  the  fabrication.  Selection  of  a  gate  metal  must  be  an  experimental 
procedure.  There  is  very  little  information  available  on  Schottky  contacts 
on  the  quaternary  alloys,  and  more  research  in  this  area  is  required. 

The  epitaxial  layer  thickness  and  impurity  doping  density  are  selected 

based  on  experience  with  GaAs  MESFETs.  The  best  devices  have  (N.)  X  thick- 

a 

12  -2 

ness  (a)  products  of  about  2.4  x  10  cm  .  Highly  doped  active  layers  are 


preferred  because  of  the  ease  of  forming  drain  and  source  ohmic  contacts 
and  because  of  the  reduction  in  the  semiconductor  resistivity.  A  reduction 
in  resistivity  produces  lower  noise  figures.  A  limit  to  the  doping  concen¬ 
tration  is  dependent  on  technological  considerations  and  the  reduction  in 
the  gate-to-channel  breakdown  voltage.  Practical  doping . concentrations  are 

in  the  range  of  5  *  10^  to  4  *  10^  cm  ®  in  GaAs  MESFET  epitaxial  layers. 

17  -3 

The  X  band  device  is  designed  with  *  1  *  10  cm  and  a  *  0.25  urn;  the 

17  -3 

Ka  band  device  with  N,  =  2.4  x  10  cm  and  a  *  0.1  pm. 

a 

The  center-fed  gate  geometry  (Figures  15  and  18)  was  selected  in  order 
to  maximize  the  RF  power  handling  capability  and  the  gain  while  maintaining 
a  low  input  impedance.  The  source  pad  area  was  maximized  to  provide  a  low 
inductive  path  to  RF  ground..  Any  source  inductance  will  provide  feedback 
that  will  increase  the  device  noise  figure  and  decrease  stability. 

The  bonding  pads  for  the  source,  drain,  and  gate  are  large  enough  to 
accommodate  0.7-mil  gold  bonding  wires  located  far  enough  away  from  the 
gate  area  so  that  the  active  region  will  not  be  damaged  by  thermocompression 
bonding.  The  source  and  drain  ohmic  contacts  are  to  be  formed  from  a 
AuGe/Ni/Au  metal  system.  The  bonding  pads  will  be  formed  by  a  gold  film 
overlay. 

The  overall  thickenss  of  the  MESFET  chips  are  selected  to  be  150  pm 
and  100  pm  for  the  X  and  Xa  band  devices,  respectively.  Mechanical  hand¬ 
ling  and  parasitic  capacitance  considerations  are  factors  in  the  selection 
of  these  particular  values.  Semi-insulating  In?  is  used  as  the  substrate 
material,  and  the  bulk  resistivity  should  have  a  minimum  value  of  about 
10®  ft* cm. 
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III.  EQUIVALENT  CIRCUIT  AND  ELEMENT  VALUE  CALCULATIONS 


An  equivalent  circuit  that  can  be  used  to  predict  the  small-signal 

operation  of  the  MESFET  devices  is  shown  in  Figure  21.  This  circuit  is 

based  on  the  model  of  Pucel,  Haus,  and  Statz.  The  circuit  element  values 

are  listed  in  Table  11  and  were  calculated  in  the  following  manner. 

1.  The  gate  resistance,  R  ,  is  due  to  the  geometry  and  resistivity 

S 

of  (he  gate  metallization.  It  has  been  shown  [35]  that  the  gate  RF 
resistance  of  short  gate  devices  is  one- third  the  dc  resistance.  There¬ 
fore, 


R 


g  3Nt  L 
g  g 


where  p  »  the  resistivity  of  the  gate  metallization  (for  Au,  p  = 

g 

3(10)"6  fl-cm) 

Z  *  gate  width 

t  »  gate  thickenss 
g 

Lg  *  gate  length 

N  »  number  of  gate  fingers  (N  =  2) . 

2.  The  source  and  drain  resistances,  R  and  R, ,  are  the  sums  of  all 

s  d’ 

the  resistance  contributions  from  the  bonding  pads  to  the  edges  of  the 
gate  depletion  region.  These  resistances  include  contributions  from  the 
contact  resistance  between  the  metal  and  the  epitaxial  layer  and  the  sheet 
resistance  of  the  semiconductor  material  between  the  edge  of  the  ohmic 
contacts  and  the  gate.  The  resistance  values  can  be  calculated  from  the 
equations 


R 


L 

_gs 

Z 


+ 


R 

c 

Z 
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Table  11.  Equivalent  Circuit  Element  Values 


• 

X  Band  Device 

Ka  Band  Device 

Rg<n) 

7.5 

12.5 

Rd  (0) 

4.80 

9.21 

R  (0) 
s 

3.13 

9.21 

Cgs  (PF) 

0.439 

0.0915 

°8d  <"r> 

0.036 

0.0086 

Cds  <PF> 

0.061 

0.021 

Rds  <s» 

762 

1142 

gm  (mmho) 

66 

55 

and 


L.  R 

R  =  p  — ^  +  — 

d  s  Z  Z 


where  p  is  the  semiconductor  sheet  resistivity  (listed  in  Table  9) ,  L 

s  gs 

and  are  the  spacings  between  the  gate  and  the  source  and  drain  elec¬ 
trodes,  and  R  is  the  contact  resistance, 
c 

The  contact  resistance,  Rc,  is  not  known  for  the  quaternary  layers 
and  is  estimated  to  be  approximately  0.05  0*cm.  Threr  could  be  consider¬ 
able  error  in  this  number,  and  the  exact  value  needs  to  be  experimentally 
established. 

3.  The  output  resistance,  R^,  is  the  resistance  between  the  drain 
and  source  with  the  channel  operating  in  velocity  saturation.  It  can  be 
calculated  from  the  expression 
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"Ivds 


4E  e_e  v  Z 
s  R  o  s 


where  V 


ds 

E_ 


'R 


'R 


drain-to-source  voltage 
saturation  field 

semiconductor  relative  permittivity 
carrier  velocity  at  the  bias  field 
drain  bias  current 
drain  current  with  no  gate  bias. 


It  has  been  shown  [36]  that  for  optimum  low-noise  operation  the 
device  should  be  biased  so  that 


The  drain-to-source  voltage  is  selected  so  that  the  bias  field  in  the 
channel  is 

E  »  8  kV/cm  , 
s 


and  this  results  in  a  carrier  velocity  of 

v  •  2  x  107  cm/sec  . 
s 

4.  The  drain-to-source  capacitance,  C^g,  is  primarily  a  function  of 
the  geometry  between  the  source  and  drain  electrodes.  It  can  be  calculated 
from  the  expression 
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K(kds) 


where 


r  a» + ^  *  Lsd>Lsd  i 

'  K  +  Lsd)ftd  +  Lsd>J 


and  Lg,  are  the  effective  source  and  drain  electrode  lengths  and  Lg<j  is 
the  source-to-drain  spacing.  The  functions  K  are  complete  elliptic  inte¬ 
grals  of  the  first  kind. 

5.  The  feedback  capacitance,  C  j,  can  be  calculated  in  the  same 
manner.  That  is, 

(cR+  Dc^Kkl  -  k^)1/2] 

Cgd  '  K(kgd) 


where 


and  is  the  drain- to-gate  spacing. 

6.  The  gate- to- source  capacitance,  C  ,  is  particularly  important 

© 

since  it  significantly  affects  the  device  frequency  response.  It  can  be 


calculated  from  the  expression 


C  »  e_e  Z 
gs  R  o 


+ H 


where  a  is  the  active  layer  thickness.  This  expression  shows  the  advantage 
of  a  short  gate  length  and  of  a  low  1^/^  ratio  for  minimizing  C^g.  The 
last  term  in  the  expression  describes  the  parasitic  capacitance  that 
results  from  the  curvature  of  the  depletion  region  at  the  edge  of  the  gate. 

7.  The  gain  of  the  defice  is  dependent  upon  the  MESFET  transconduc¬ 
tance,  g^.  This  parameter  can  be  calculated  with  the  equation 


2e_e  v  Zf  . 

R  o  s _ 1 

!m  a  1  -  I./I 

L  d  s 


As  can  be  seen,  a  high  v  ,  large  Z,  and  small  a  are  desirable  to  produce  a 

s 

large  transconductance.  The  bias  conditions  also  have  a  significant  effect 
upon  this  parameter,  and  for  optimum  gain,  Id/Ig  should  be  maximized.  Such 
a  bias  condition,  however,  results  in  an  increase  in  the  device  noise  figure, 
and  the  requirements  for  maximum  gain  and  low-noise  operation  are  in  con¬ 
flict. 

8.  The  inductances  L  ,  L  ,  and  L ,  represent  the  bond  wires  that  must 

x  g  d 

be  used  to  interface  the  device  with  the  RF  circuit.  If  0.7-mil  gold  wire 
is  used,  the  inductances  can  be  calculated  from  the  relationship 


50  mil  3  1  nh  . 


Parallel  wires  would,  of  course,  reduce  the  parasitic  inductance. 


IV.  FIGURES-OF-MERIT  AND  BIAS  CONDITIONS 


The  anticipated  performance  of  these  devices  are  indicated  by 
calculation  of  certain  figures-of -merit .  The  calculated  numbers  are 
listed  in  Table  12  and  were  determined  in  the  following  manner. 


(a)  Cutoff  frequency  is  the  frequency  at  which  the  forward  current 
gain  is  unity.  It  can  be  thought  of  as  the  gain-bandwidth 
product  and  can  be  expressed  as 

.  «m 

fT  "  2irC 

gs 

(b)  The  maximum  frequency  of  oscillation  is  the  frequency  at  which 
the  device  matched  power  gain  has  decreased  to  unity  and  is 
expressed  as 


max 


2^+f^ 


where 


R  +R 
g  s 

1  "  R. 

ds 


and 


T  -  2ttR  C . 

g  dg 


(c)  Mason's  unilateral  gain  is  approximately 

i2 


°  ■  w 
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It  is,  of  course,  desirable  to  maximize  all  of  these  figures-of- 
merit  for  optimum  RF  performance. 

The  ^ESFET  bias  conditions  were  also  calculated  and  are  listed 
in  Table  13.  The  calculated  values  were  determined  in  the  following 
manner. 

(d)  The  drain  current  with  zero  gate  bias,  I  ,  is  given  by  the 

s 

expression 

I  *  a  qviN.ZE 
s  n  d  s 

where  pm  is  the  low-field  mobility. 

(e)  The  dc  bias  current  was  determined  from  the  low-noise  bias 
condition  that 

I,  «  0.15  I  . 
d  s 

(f)  The  drain-to-source  bias  voltage  and  the  gate  voltage  are  to  be 

adjusted  to  provide  the  correct  bias  current.  An  estimate  for  V, 

as 

is  obtained  by  requiring  the  channel  field  to  be  sufficient  for 
velocity  saturation.  The  channel  field  was  selected  to  be  8  kV/cm. 

(g)  The  gate  pinch-off  voltage  can  be  determined  from  the  expression 

qNda2 

^pinch-off  ~  2e_e 
R  o 
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Table  12.  Calculated  Figures-of-Merit 


X-Band  Device  K  -Band  Device 

_ _ a  _ 

Cutoff  Frequency,  f^,  23.9  GHz  95.7  GHz 

Max.  frequency  of  51.2  GHz  165  GHz 

Oscillation,  f 

max 

Unilateral  Gain,  U  12.6  db  @  12  GH  12.3  db  @  40  GHz 


Table  13.  Device  dc  Parameters 


X-Band  Device 

K  -Band  Device 
a 

427 

136 

lia  <-t> 

64 

20.5 

Vds  (v) 

4 

2 

7pinch-off  ^ 

4.28 

1.64 
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V.  SUMMARY 


Designs  for  MESFETs  capable  of  X-  and  Ka-band  operation  have 
been  presented.  The  X-band  device  utilizes  standard  one-micron 
gate  length  technology.  An  optimized  Ka-band  MESFET,  however,  will 
require  state-of-the-art  fabrication  technology.  The  devices  are 
intended  to  take  advantage  of  the  favorable  material  parameters  of  the 
quaternary  alloy  Ga^  27In0  73F0  4As0  6  mater^a^  system  lattice  matched 
to  semi-insulating  InP  substrates .  Such  devices  should  provide  improved 
performance  over  comparable  GaAs  devices. 

An  equivalent  circuit  based  upon  the  model  of  Pucel,  Haus  and 
Statz  has  been  derived  for  the  devices.  The  equivalent  circuit  can  be 
used  for  predictions  of  device  performance.  Much  of  the  MESFET  design 
was  based  upon  the  optimization  of  the  model  circuit  elements. 

An  indication  of  device  performance  is  given  by  the  calculated 
figures-of -merit.  These  calculations  indicate  that  the  MESFET 
devices  have  considerable  potential  as  low-noise,  high-gain  microwave 
and  millimeter  wave  amplifiers. 
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